Experimental determination of the Berry phase in a superconducting charge pump 
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We present the first measurements of the Berry phase in a superconducting Cooper pair pump. 
A fixed amount of Berry phase is accumulated to the quantum mechanical ground state in each 
adiabatic pumping cycle, which is determined by measuring the charge passing through the device. 
The dynamic and geometric phases are identified and measured quantitatively from their different 
response when pumping in opposite directions. Our observations, in particular the dependencies of 
the dynamic and geometric effects on the superconducting phase bias across the pump, agree with 
the basic theoretical model of coherent Cooper pair pumping. 



Geometric phases arise from adiabatic cyclic evolution 
in classical and quantum physics In contrast to dy- 
namic effects, geometric phases depend only on the geom- 
etry of the cycle traversed by the state of the system. In 
parallel transport, for example, a vector is moved along 
a path without changing its direction in a local coordi- 
nate system For a closed path in a flat Euclidean 
space, the vector returns exactly to its initial state with 
respect to a global coordinate system. In a curved space 
however, the direction of the vector can change as shown 
in Fig. [Ha), where the vector moves on a surface of a 
sphere along a loop enclosing a solid angle fl. In this 
case, the angle 9 between the initial and the final state 
of the vector equals f2 which depends only on the cho- 
sen path, and hence is regarded as a geometric phase. In 
our universe, geometric phases have been employed for 
example to measure the curvature of space due to grav- 
itation, and hence to test Einstein's theory of general 
relativity This geodetic effect is intended to be mea- 
sured as changes in the rotation axes of gyroscopes inside 
a satellite orbiting Earth with the stringent accuracy of 
10~ 4 in the project Gravity Probe B Q. 

We consider adiabatic and cyclic temporal evolution in 
the ground state of a quantum-mechanical system. The 
state of any pure quantum system can be described by a 
complex valued wave function. Thus the simplest ge- 
ometric phase accumulated in a cycle, i.e., the Berry 
phase [4j, is a phase shift of the complex number mul- 
tiplying the wave function. As such, the absolute phase 
of the wave function does not have a physical meaning, 
and hence is unobservable. Thus the measurements of 
the Berry phase typically rely on the interference of two 
states which have undergone a different phase shift [j], [|| . 
This technique is also employed in the proposal to mea- 
sure the Berry phase in an asymmetric superconducting 
quantum interference device (SQUID) In phase bi- 
ased Cooper pair pumps 0] however, the accumulated 
Berry phase is related to the pumped charge d, Q , and 
hence we have a fundamentally different way to deter- 
mine it. Here, we report on the first experimental real- 
ization of phase biased Cooper pair pumping in a super- 
conducting loop. 
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FIG. 1: (a) Parallel transport of a vector along a path (red 
line) enclosing a solid angle f2. The lightest arrow shows the 
initial state of the vector and the darkest arrow the final state. 
The resulting angle between the initial and the final vectors 9 
is equal to the solid angle £1. (b) Scanning electron micro- 
graph of the island on the left and of the detector on the 
right, (c) Simplified circuit diagram of the measured sample. 
The corresponding parts in the circuit diagram and SEM- 
images are marked by colors. The independent fluxes <J?i, $2, 
and $dc are controlled by on-chip coils and the gate volt- 
age V g is related to the gate charge n g and gate capacitance C g 
by n g = C g V g /(2e). 



Recently, superconducting circuits have proven to 
be suitable for coherent manipulations of quantum 
states [H|,[l2|]) m particular, two-level quantum systems. 
The natural scalability of electric circuits makes them 
potential candidates for qubits, i.e., basic building blocks 
of the emerging quantum computer. On the other hand, 
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holonomies arising from adiabatic and cyclic evolution 
in a degenerate eigenspace 13( are unitary transforma- 
tions which in turn, can be compiled to execute quan- 
tum algorithms 1J] of practical interest. In particular, 
the holonomies related to charge transport in Cooper pair 
devices have been studied theoretically [l5[ . Thus the ob- 
servation of the Berry phase in superconducting circuits 
is an important step towards the development of holo- 
nomic quantum computation To date, holonomic 

quantum computation has only been demonstrated us- 
ing liquid state nuclear magnetic resonance [ijj (NMR), 
the scalability of which is limited to about ten qubits. 

Single-charge pumping is based on transporting a con- 
trolled number n of carriers with quantized charge e* in a 
cycle repeated at frequency /. This principle yields ide- 
ally the pumped charge Q p = ne* , and hence the pumped 
current 



I p = ne*f. 



(1) 



The carrier charge e* is e for single electrons and 2e for 
Cooper pairs. In charge pumps described by a phase 
coherent order parameter field, the phase difference of the 
field across the device ip may play a significant role ■ In 
fact for a constant ip, adiabatic charge pumping gives rise 
to the Berry phase [1, @] . Interestingly, the accumulated 
Berry phase, Ob, in a cycle is related to the pumped 
charge, Q p , by 



Qp = -e*dO B /dip. 



(2) 



In an ideal pumping cycle corresponding to Eq. (p}, the 
accumulated Berry phase equals — nip. However, this 
kind of pumping does not yield definite fingerprints on 
the relation between the Berry phase and the pumped 
charge. Therefore, it is important to reveal the phase 
coherent nature of pumping from its dependence on ip. 
Allowing a non- vanishing average persistent current, i.e., 
leakage, through the pump during the cycle, the pumped 
charge becomes phase dependent. In the regime where 
charge states of the pump are approximately the eigen- 
states of energy, one obtains in the two-charge-stat e ap - 
proximation for a cycle described in Fig. [2ja) @, @, llq 



Q p w e*n(l — 5 cos ip), 



(3) 



where (5 <gC 1 is proportional to the leakage of the pump. 
Due to its geometric origin, the pumped charge is inde- 
pendent of the pumping frequency. The charge leaked 
through the pump during a cycle, i.e., the dynamic part 
of the transferred charge is obtained as 
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where T = 1// is the period of the pumping cycle, / cma 
is the maximum critical current in the pumping cycle, 
and P is a constant specific to the particular pump and 
to the control parameter cycle. To justify the validity of 



the above model in our measurements, we compare the 
measured pumped current I p = Q p f and the dynamic 
current Id = Qdf with Eqs. ([3]) and ([4|). Here, Id is the 
supercurrent in the ground state of the system averaged 
over the pumping path and thus it does not depend on / 
or n, nor on the direction of traversing the pumping path. 
On the contrary, I p is proportional to both / and n in the 
adiabatic evolution, and its sign changes on reversing the 
path. Agreement with the theory allows us to determine 
the accumulated Berry phase from Eqs. (J2J) and |3j) as 
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—n(ip — 5 shop). 



(5) 



Figure QJb) and (c) shows the Cooper pair pump, the 
sluice in a 800 am 2 superconducting loop with a de- 
tector junction. In the ideal pumping cycle shown in 
Fig.[2ja), an integer number of n™ ax excess electron pairs 
is first attracted to the island through one SQUID and 
then repelled from the island through the other one us- 
ing the gate voltage and the tunable critical currents of 
the SQUIDs in analogy with a piston pump. Hence the 
sluice generates ideally an average pumped current given 
by Eq. |T]) with n = n™ ax The details of the working 
principle of the sluice can be found in Refs. HGH. 

The Josephson junctions denoted by black crosses in 
the circuit diagram of Fig. QJc) consist of A10 x tunnel 
barriers fabricated by standard electron beam lithogra- 
phy and two-angle evaporation into an all-aluminum de- 
vice on oxidized silicon wafer. The sluice part of our 
sample is identical to the one used in Ref. [10( except for 
up to 10% smaller junction size and stronger oxidation. 
The charging energy of the island is difficult to measure 
in the presence of the detector junction, and hence it is 
estimated based on the sample used in Ref. [nj to be 
2 K x ks ■ Using the Ambegaokar-Baratoff formula and 
the IV characteristics of the sample, the maximum crit- 
ical currents of the SQUIDs and the detector were esti- 
mated to be 30 nA each and 70 nA, respectively. Due to 
the parallel structure of the junctions in the circuit, there 
is an uncertainty of up to 20% in the parameter estima- 
tion. The plasma frequency of the detector is roughly 
20 GHz. 

We install a Josephson junction working as a thresh- 
old current detector in parallel with the sluice forming 
a superconducting loop as proposed in Ref. [8J]. Fig- 
ure[ljc) also shows our measurement scheme, in which we 
feed current pulses through the circuit and monitor volt- 
age across it. A voltage pulse is observed if the system 
switches into the normal state in response to the cur- 
rent pulse. The repetition rate of the measurement was 
adjusted low enough for the switching events to be uncor- 
rected. Such a detection method of circulating current 
was realized, e.g., in the measurements of the supercon- 
ducting qubit Quantronium To assure feasible op- 
eration, the critical current of the detector is chosen to 
be much higher than any instantaneous critical current 
of the sluice. The probability for the system to switch to 
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FIG. 2: (a) Ideal pumping cycle (solid line) in the control 
parameter space of the sluice. The tunable critical currents 
of the SQUIDs in the sluice 7 c i = 7 c i (<&i) and J C 2 = I C 2(&2) 
are modulated together with the gate charge at the pumping 
frequency / = 14 MHz for all the data in this paper. In this 
cycle, the two critical currents are analogous to valves of the 
charge flow, and the gate charge to a piston. The dashed 
line illustrates a pumping cycle with finite J crcs = <5/ cmax /2 
introducing finite dynamic current, (b) Pumped current as 
a function of the maximum gate charge n™ ax . The blue line 
shows the theoretical value obtained from Eq. ([1]) which is 
expected to be valid in the case of adiabatic pumping for 
small pumping amplitudes n™ ax . There are no fitted param- 
eters since the conversion constant of the gate voltage to gate 
charge is obtained from DC measurements. 



the normal state depends strongly on the height of the 
external current pulse / pu iso- The switching current /50 of 
the system is defined to be the point where the probabil- 
ity is 50%. We define backward pumping as the direction 
for which the pumped current adds to the applied 400 [is 
current pulse, and hence the current through the detector 
is given by I d ot = ipuise + I p — Id- For forward pumping 
obtained by traversing the pumping cycle in the opposite 
direction, the pumped current compensates part of the 
external pulse: idet = -fpuisc — Ip — Id- Thus we measure 
shifts in /50 which correspond to twice the pumped cur- 
rent. The average of the two switching currents equals 
to the dynamic current plus a constant, namely, T^et at 
the switching point. Note that we restrict our studies to 
take into account only the dc component of the pumped 
current at the detector. The effect of the ac component 
not filtered by the circuit is left for future research. 

Figure H[b) displays the measured pumped current as 
a function of the gate amplitude for our most ideal pump- 
ing cycle. For low enough gate amplitudes, the data 
shows a nearly linear dependence in good correspondence 
with the theoretical behaviour given in Eq. (JTJ) . For high 
gate amplitudes, the adiabaticity of the pump breaks 
down and deviation from Eq. |T]) is observed. These re- 
sults demonstrate the first observation of Cooper pair 
pumping in closed superconducting circuits [8|, [19( . How- 
ever, the data in Fig.^b) does not prove that the current 
arises from coherent quantum dynamics. 

The phase difference ip of the superconductor order pa- 



rameter across the sluice has to be a classical parameter 
for Eq. f2| to hold. This is satisfied due to the detector 
junction which protects the sluice from voltage fluctua- 
tions, and hence phase biases it. For a detector junction 
with large enough critical current, the phase difference 4>d 
across it is obtained from arcsin(idet /Ic) at the switch- 
ing point where / pu isc = ho- Thus we can control ip by 
adjusting magnetic flux $dc through the loop and using 
the fundamental phase relation of a superconducting loop 
(p-(f> d = 27t$dc/*o, where $ = h/(2e) = 2.07 fWb is 
the flux quantum. Due to proper magnetic shielding, the 
flux offset in $dc was negligible. 

The evidence of the phase coherence in our measure- 
ment is shown in Fig. [3j Figure [3ja) presents the varia- 
tion of the switching current of the system with respect 
to the phase difference across the sluice controlled by the 
eternal flux <I>dc as described above, and Fig. [3jb) dis- 
plays the pumped current from the same measurement. 
Clear sinusoidal modulation is observed in both curve 
sets in agreement with Eqs. © and (|3]), implying that 
the sluice is coherent and phase biased. Note that the 
minimum of the pumped current corresponds quite ac- 
curately to the point of vanishing phase difference across 
the sluice and the dynamic current is phase shifted by 
almost 7r/2 radians as in Eqs. © and Q. The pumped 
number of Cooper pairs n and the relative modulation 
amplitude 5 are determined from the data and the cor- 
responding curves for the accumulated Berry phases are 
shown in Fig. (He) according to Eq. ([5]). For the largest 
values of 6, second order corrections to Eq. ([5]) may mod- 
ify the estimated value for the Berry phase. 

To further test our scheme, we measured the depen- 
dence of the modulation amplitudes on the gate ampli- 
tude and residual critical current 7 crcs = <5/ C max/2, i.e., 
leakage, shown in Figs.[3{d) and[3{e). The residual crit- 
ical current was introduced by an offset in the control 
fluxes of the SQUIDs from the ideal pumping cycle. Fig- 
ure[3jf) shows S as a function of the offset. As predicted, 
the modulation amplitude of both, the pumped and the 
dynamic current, increase with the residual critical cur- 
rent, but only the pumped current depends on the maxi- 
mum gate charge n™ ax , in agreement with the theory. We 
note that this way of measuring the Berry phase is funda- 
mentally different from the conventional method utilizing 
interference with excited states. 

Our observations pave the way for further experiments 
on Cooper pair pumping in closed circuits (l9(, on the 
quantum standard of electric current, for applications 
of geometric phases in holonomic quantum computa- 
tion . and test the fundamental implications of the 
quantum theory in an order parameter describing co- 
herent dynamics of a macroscopic number of condensed 
Cooper pairs. Our measurements are in agreement with 
a theoretical model neglecting decoherence. [? ] . On 
the other hand, the effects of dephasing and dissipation 
in superconducting qubits are harmful in the manipu- 
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FIG. 3: (a) Typical switching currents of the system averaged 
over forward and backward pumping directions as a function 
of the phase ip across the sluice for several residuals of the 
critical current Teres, see Fig. [2ja). Here, n™ ax = 16. (b) 
Same as panel (a) but for the pumped current. To obtain 
the amount of accumulated Berry phase, we make a fit to 
the data corresponding to Eq. © with 5 and n as fitting 
parameters, the result of which is shown in panel (c). (d) The 
variation amplitude of the switching current as in panel (a) 
with respect to the phase across the sluice as a function of the 
offset in the RF fluxes introducing residual critical current for 
n max = 4 ( b i ue ^ 10 (green), 16 (red), 22 (cyan), 28 (magenta), 
35 (yellow), and 41 (gray). These values of the maximum 
gate charge are in the linear adiabatic regime as shown in 
Fig-Efb). (e) The variation amplitude of the pumped current 
as in panel (b) from the same data set as in panel (d). (f) 
Modulation amplitude in the pumped current divided by the 
average pumped current (shown in the inset), 8, as a function 
of the flux off-set. The dashed line shows a linear fit to the 
data for eight smallest off-sets, and hence yields the linear 
dependence of the leakage parameter 8 on the flux off-set. The 
dashed line in panel (d) shows a corresponding theoretical line 
from Eq. © using f3 = 0.0245 and L„ = 30 nA. The dashed 
lines in panel (e) show the linear dependence of 8 multiplied 
by the average pumped current (see the inset) corresponding 
to the specific gate amplitude and the smallest flux off-set. 



lation of superconducting quantum systems, see, e.g. 



Refs. [2l| . Thus our results support the robustness of geo- 
metric phases against decoherence and fluctuations [22j . 
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— Note added. During the peer review of this work 
Leek ct al. reported on the observation of the Berry phase 
in a superconducting qubit 23]. 
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